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Abstract 

& : 

We present a path-integral formulation of 't Hooft's derivation of quantum from 
classical physics. The crucial ingredient of this formulation is Gozzi et al. super- 
symmetric path integral of classical mechanics. We quantize explicitly two simple 
classical systems: the planar mathematical pendulum and the Rossler dynamical 
system. 
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In recent decades, various classical, i.e., deterministic approaches to quantum theory have been 
proposed. Examples are Bohmian mechanics [1], and the stochastic quantization procedures of 
Nelson [2], Guerra and Ruggiero [3], and Parisi and Wu [4,5]. Such approaches are rinding in- 
creasing interest in the physics community. This might be partially ascribed to the fact that 
such alternative formulations help in explaining some quantum phenomena that cannot be easily 
explained with the usual formalisms. Examples are multiple tunneling [6], critical phenomena at 
zero temperature [7] , mesoscopic physics and quantum Brownian oscillators [8] , and quantum- field- 
theoretical regularization procedures which manifestly preserve all symmetries of the bare theory 
such as gauge symmetry, chiral symmetry, and supersymmetry [9]. They allow one to quantize 
gauge fields, both Abelian and non-Abelian, without gauge fixing and the ensuing cumbersome 
Faddeev-Popov ghosts [10], etc.. 

The primary objective of a reformulation of quantum theory in the language of classical, i.e., 
deterministic theory is basically twofold. On the formal side, it is hoped that this will help in 
attacking quantum-mechanical problems from a different direction using hopefully more efficient 
mathematical techniques than the conventional ones. Such techniques may be based on stochastic 
calculus, supersymmetry, or various new numerical approaches (see, e.g., Refs. [5,11] and citations 
therein). On the conceptual side, deterministic scenarios are hoped to shed new light on some old 
problems of quantum mechanics, such as the origin of the superposition rule for amplitudes and 
the theory of quantum measurement. It may lead to new ways of quantizing chaotic dynamical 
systems, and ultimately a long-awaited consistent theory of quantum gravity. There is, however, 
a price to be paid for this; such theories must have a built-in nonlocality to escape problems with 
Bell's inequalities. Nonlocality may be incorporated in numerous ways — the Bohm-Hiley quantum 
potential [1,12], Nelson's osmotic potential [2], or Parisi and Wu's fifth-time parameter [4,5]. 

Another deterministic access to quantum-mechanical systems was recently proposed by 't Hooft 
[13,14] with subsequent applications in Refs. [15-21]. It is motivated by black-hole thermodynamics 
(and particularly by the so-called holographic principle [22,23]), and hinges on the concept of 
information loss. This and certain accompanying non-trivial geometric phases are able to explain 
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the observed non-locality in quantum mechanics. The original formulation has appeared in two 
versions: one involving a discrete time axis [16], the second continuous times [14]. The goal of 
this paper is to discuss further and gain more understanding of the latter model. The reader 
interested in the discrete-time model may find some practical applications in Refs. [24,25]. It is 
not our purpose to dwell on the conceptual foundations of 'tHooft's proposal. Our aim is to set 
up a possible useful alternative formulation of 'tHooft's model and quantization scheme that is 
based on path integrals [11]. It makes use of Gozzi et al. path- integral formulation of classical 
mechanics [26,27] which appears to be a natural mathematical framework for such a discussion. 
The condition of the information loss, which is basically a first-class subsidiary constraint, can 
then be incorporated into path integrals by standard techniques. Although 'tHooft's procedure 
differs in its basic rationale from stochastic quantization approaches, we show that they share a 
common key feature, which is a hidden BRST invariance, related to the so-called Nicolai map [28]. 
To be specific, we shall apply our formulation to two classical systems: a planar mathematical 
pendulum and the simplest deterministic chaotic system — the Rossler attractor. Suitable choices 
of the "loss of information" condition then allow us to identify the emergent quantum systems 
with a free particle, a quantum harmonic oscillator, and a free particle weakly coupled to Duffing's 
oscillator. 

Our paper is organized as follows. In Section II we quantize 'tHooft's Hamiltonian system 
by expressing it in terms of a path integral which is singular due to the presence of second- 
class primary constraints. The singularity is removed with the help of the Faddeev-Senjanovic 
prescription [29,30]. It is then shown that the fluctuating system produces a classical partition 
function. In Section III we briefly review Gozzi et al. path-integral formulation of classical 
mechanics in configuration space. The corresponding phase-space formulation is more involved 
and will not be considered here. By imposing the condition of a vanishing ghost sector, which is 
characteristic for the underlying deterministic system, we find that the most general Hamiltonian 
system compatible with such a condition is the one proposed by 'tHooft. In Section IV we 
introduce 'tHooft's constraint which expresses the property of information loss. This condition 
not only explicitly breaks the BRST symmetry but, when coupled with the Dirac-Bergmann 
algorithm, it also allows us to recast the classical generating functional into a form representing 
a proper quantum-mechanical partition function. Section V is devoted to application of our 
formalism to practical examples. We conclude with Section VI. For the reader's convenience the 
paper is supplemented with four appendixes which clarify some finer mathematical points needed 
in the paper. 



II. QUANTIZATION OF 'THOOFT'S MODEL 

Consider the class of systems described by Hamiltonians of the form 

N 

H = Y,Pafa{ci) ■ (1) 

o=l 

Such systems emerge in diverse physical situations, for example, Fermi fields, chiral oscillators [20], 
and noncommutative magnctohydrodynamics [31]. The relevant example in the present context 
is the use of (1) by 'tHooft to formulate his deterministic quatization proposal [13]. 

An immediate problem with the above Hamiltonian is its unboundedness from below. This is 
due to the absence of a leading kinetic term quadratic in the momenta p 2 a /2M, and we shall dwell 
more on this point in Section IV. The equations of motion following from Eq.(l) are 

qa = fa(q), Pa = -Va—R • (2) 

Note that the equation for q a is autonomous, i.e., it is decoupled from the conjugate momenta 
p a - The absence of a quadratic term makes it impossible to find a Lagrangian via a Legendrc 
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transformation. This is because the system is singular — its Hess matrix H ah = d 2 H / dp a dp 
vanishes. 

A Lagrangian yielding the equations of motion (2) can nevertheless be found, but at the expense 
of doubling the configuration space by introducing additional auxiliary variables q a (a — 1, . . . , N). 
This extended Lagrangian has the form 

N 

L = ^2[q a q a - <?a/a(q)] (3) 

a=l 

and it allows us to define canonically conjugate momenta in the usual way: p a = dL/dq a , p a = 
dL/dq a . A Legendre transformation produces the Hamiltonian 

N N 

H(p a ,qa,Pa,q~a) = ^PaQa+PaQa ~ L = ^<L/a(q) • (4) 
a— 1 a— 1 

The rank of the Hess matrix is zero which gives rise to 2N primary constraints, which can be 
chosen as: 

4> a i=Pa-q a « 0, <f%=p a W 0. (5) 

The use of the symbol w instead of = is due to Dirac [32] and it has a special meaning: two 
quantities related by this symbol are equal after all constraints have been enforced. The system 
has no secondary constraints (see Appendix A) . The matrix formed by the Poisson brackets of the 
primary constraints, 

{#(*),$(*)} = Sab, (6) 

has a nonzero determinant, implying that all constraints are of the second class. Note that on the 
constraint manifold the canonical Hamiltonian (4) coincides with 'tHooft's Hamiltonian (1). 

To quantize 'tHooft's system we utilize the general Faddeev-Senjanovic path integral for- 
mula [29,30] for time evolution amplitudes 1 

(q2,t 2 |qi,ti) = AT J VpVq ^J\det |{&, ^}||| exp j 1 jf ' dt [pq - H (q, p)] J . (7) 

Using the shorthand notation fa — fa], fa], fa\,fa\, ■ ■ ■ , 4>i , faf (* = 1j • • • i 2^V), Eq.(7) implies in 
our case that 

(q2,*2|qi,*i) = Af J VpVqVpVq S[p - q] <5[p] expj^ dt [pq + pq H (q, q, p, p) 

/■q(*2)=q2 r j rt 2 

= Af / 2?qX>q exp - / L(q,q,q,q) 

>/q(ti)=qi L ft Jti 

/•q(*2)=q2 

= aa/ pq n%-/.(q)l. (8) 

Jq(ti)=qi „ 

where 5[f] = 5(f(i)) is the functional version of Dirac's 5-function. This result shows that 
quantization of the system described by the Hamiltonian (1) retains its deterministic character. 
The paths are squeezed onto the classical trajectories determined by the differential equations 



1 Other path-integral representations of systems with second-class constrains such as that of Fradkin and 
Fradkina [33] would lead to the same result (8) . 



3 



ia = /a(q)- The time evolution amplitude (8) contains a sum over only the classical trajectories 
- there are no quantum fluctuations driving the system away from the classical paths, which is 
precisely what we expect from a deterministic dynamics. 

The amplitude (8) can be brought to a more intuitive form by utilizing the identity 

5[f(q)-q] = %-q c i] (detM)- 1 , (9) 

where M is a functional matrix formed by the second derivatives of the action -4[q, q] = 
JdtL(q, q,q, q) : 



M ab (t,t r ) 



S 2 A 



Sq a (t) 6q b (t>) 



(10) 



The Morse index theorem then ensures that for sufficiently short time intervals t 2 — t\ (before the 
system reaches its first focal point), the classical solution with the initial condition q(ti) = qi is 
unique. Note, however, that because of the first-order character of the equations of motion we are 
dealing with a Cauchy problem, which may happen to possess no classical trajectory satisfying 
the two Dirichlet boundary conditions q(t\) = qi, q(^) = q2- If a trajectory exists, Eq. (8) can 
be brought to the form 



/■q(*2)=q2 

(q2,*2|qi,ii) = M / 2?q (5 [q - q c[ ] , 

Jq(tl)=qi 



(11) 



where M = A/*/(detM). We close this section by observing that detM can be recast into more 
expedient form. To do this we formally write 



detM = det 



cxp 



dtS b n + 



Tr In a 



6 , df a (q(t)) 



dq b (t) 



5(t-t') 



= exp 



Trln 



6 , 0/a(q(*)) 



dtS b n + 



dq b (t) 



S(t - t') 



S b a S(t - t') + G(t - t') 



J/.(qW) 



dq b (t') 



exp [Tr(ln<9t)] cxp 



Trln 



5 b J(t f) + G(t t>) 



dq b (t') 



(12) 



Here G(t — t') is the Green's function satisfying the equation 

d t G(t-t') = S(t-t'). 

Choosing G(t — t') = 9{t — t'), and noting that the first factor in Eq.(12) is an irrelevant constant 
that can be assimilated into TV we have 



detM = exp 



= exp 



Trln 



5 b J(t f ) + G(t f ) 



,%(q(t')) 



dq b (t>) 



cxp 



Tr 



dq b (t) 



1 f t2 

- jf dt V q f (q) 



(13) 



In deriving Eq.(13) we have used the fact that due to the product of the ^-function in the expansion 
of the logarithm, all terms vanish but the first one. In evaluating the generalized function 9(x) at 
the origin we have used the only consistent midpoint rule [11]: 9(0) = 1/2. Using the identity 



exp 



1 f t2 



dt V q f (q) 



q=q c i 



J Vq S [q - q cl ] cxp ~\ ^ dt V q q 



(14) 



we can finally write the amplitude of transition in a suggestive form 



(q2»*2|qi,*i) = 

= TV 



q(*2)=q2 

q(*i)=qi 
q(*2)=q2 

q(*i)=qi 



VqVq S[q- q c i]<S[q - q c i] exp 



2 k 



dt V q q 



VqVq S[q- q c \}S[q - q cl ] 



I det K(t 2 ) 
dctif(ti) 



(15) 



Here K{t) is the fundamental matrix of the solutions of the system 



la 



-<lb 



dq a 



(16) 



det K(t) is then the corresponding Wronskian. Note that in the particular case when V q f (q) = 0, 
i.e., when the phase flow preserves the volume of any domain in the configuration space, the 
exponential in Eq.(15) can be dropped. 2 Because the exponent depends only on the end points 
of q variable it can be removed by performing the trace over q. As a result we can cast the 
quantum-mechanical partition function (or generating functional) Z into the form 



7 



Z = M \ VqDq <5[q - q c iMq - q cl ] exp 



t 2 



[J(t)q(t)+J(t)q(t)]dt 



Af / Vq a 5[q a - (q a )ci] exp 



dt J a {t)q a {t) 



(17) 



Here the doubled vector notation q a = {q, q} and J a = {J, J} was used. 



III. PATH INTEGRAL FORMULATION OF CLASSICAL MECHANICS - 
CONFIGURATION-SPACE APPROACH 

Expressions (11) and (17) formally coincide with the path- integral formulation of classical me- 
chanics in configuration space proposed by Gozzi [26] and further developed by Gozzi, Reuter, 
and Thacker [27] (see also Ref. [21] for recent applications). Let us briefly review aspects of this 
which will be needed here. Consider the path-integral representation of the generating functional 
of a quantum- mechanical system with action A[q]: 



Zqm=M J Vq e-^[ q l/ R exp J 3(t)q(t)dt 



(18) 



We assume in this context that there are no constraints that would make the measure more 
complicated as in Eq. (7). Gozzi et al. proposed to describe classical mechanics by a generating 
functional of the form (18) with an obviously modified integration measure which gives equal 
weight to all classical trajectories and zero weight to all others 



Zcm = M j Vq S[q - dci] exp J 3(t)q(t)dt 



(19) 



Although the form of the partition function (19) is not derived but postulated, we show in Ap- 
pendix B that it can be heuristically understood either as the "classical" limit of the stochastic- 
quantization partition function (c.f., Appendix BI), or as a results of the classical limit of the 
closed-time path integral for the transition probability of systems coupled to a heat bath (c.f., 
Appendix BII). This, in turn, indicates that it would be formally more correct to associate (19) 
with the probability of transition or (via the stochastic-quantization passage) with the Euclidean 
amplitude of transition [34]. Albeit (19) cannot be generally obtained from (18) by a semiclassical 
limit a la WKB (which can be recognized by the absence of a phase factor exp(z / HA(q c \)) in (19)) 
it may happen that even ordinary amplitudes of transition posses this form. This is the case, for 
instance, when the number of degrees of freedom is doubled or when one deals with closed-time- 
path formulation of thermal quantum theory. Yet, whatever is the origin or motivation for (19), 
it will be its formal structure and mathematical implications that will interest us here most. 



To proceed we note that an alternative way of writing (19) is 



2 This corresponds to the situation when there are no attractors in the configuration space T q . 
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Z C m —Ml T>q 5 



/ 



^4 
«5q 



det 



S 2 A 



5q a (t) Sq b (t') 



exp 



J 3(t)q(t)dt 



By representing the 5 functional in the usual way as a functional Fourier integral, 



SA 



Sq 



VX exp i 



ri2 

i / dt X(t) 



SA 



(20) 



(21) 



and the functional determinant as a functional integral over two real time-dependent Grassmannian 
ghost variables c a (t) and c a (t), 



det 

we obtain 



S 2 A 



5q a (t) Sq h (t') 



Zqm = J VqDXVcDc exp 
with the new action 

rt 2 



rt-2 

iS + dt J(i)q(i) 

Jtl 



(22) 



(23) 



r 2 ^ f 2 r 2 
5 [q ,e,c,A] s jf *m m -ij h *J h tfj.Wjj^j,^. (") 

Since Zqm together with the action (24) formally result from the classical limit of the stochastic- 
quantization partition function, it comes as no surprise that S exhibits BRST (and anti-BRST) 
supersymmetry. It is simple to check that S does not change under the supersymmetry transfor- 
mations 



<feRSTq = ec, <5 B rstc = 0, <S B rstc = —ieX, <5brstA = 0, 



(25) 



where e is a Grassmann-valued parameter (the corresponding anti-BRST transformations are 
related with (25) by charge conjugation). Indeed, the variations of the two terms in (24) read 



<5brsi 



->BRST 



* 2 , w x SA 



*2 



dt 



dt' X a {t) 



S 2 A 



7T Cb(t') , 



6q a (t)5q b (t') 



7T C b (t') 



5q a (t)5q b (t>) 



is dt dt' X a (t) 



/•*2 rt 2 

/ dt dt' c a 

Jtt Jt! 



Sq a (t)5q b (t') 



(26) 



7T C6(f ) 



(27) 



The second term on the RHS of (27) vanishes because the functional derivative of A is symmetric 
in c <-> 6 whereas the term c c c;, is anti-symmetric. Inserting Eqs.(26) and (27) into the action we 
clearly find <5brst<? — 0. As noted in [27], the ghost fields c and c are mandatory at the classical 
level as their role is to cut off the fluctuations perpendicular to the classical trajectories. On the 
formal side, c and c may be identified with Jacobi fields [27,35]. The corresponding BRST charges 
are related to Poincare-Cartan integral invariants [36]. 

By analogy with the stochastic quantization the path integral (23) can, of course, be rewritten 
in a compact form with the help of a superfield [26,34] 

*„(t,0,0) = qa(t)+i6c a (t)-i6c a (t)+i66X a (t), (28) 

in which 9 and 9 are anticommuting coordinates extending the configuration space of q a variable 
to a superspace. The latter is nothing but the degenerate case of supersymmetric field theory in 
d = 1 in the superspace formalism of Salam and Strathdee [37] . In terms of superspace variables 
we see that 
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J d9d9 A[&] = J dtdQdO L(q(t) + iOc(t) - i9c(t) + i96\(t)) 

= J dSdO A[q\ + J dtdSdO (iOc(t) - iSc(t) + i66X) 



/ 



dtdt'dUO 9c a (t) S *f 6c{t'). (29) 
5q a (t)Sq b (t') 



Using the standard integration rules for Grassmann variables, this becomes equal to — iS. Together 
with the identity T>$> = VqDcDcDX we may therefore express the classical partition functions 
(19) and (20) as a supersymmetric path integral with fully fluctuating paths in superspace, 



■i CM 



j V<f> cxpj-^ d6d6 A[&] (6,6) + J dtd9d9 r(t, 0, 0)*(t, 0, 0) j • (30) 



Here we have defined the supercurrent T(t, 0, 9) = 993 (t). 

It is interesting to find the most general form of an action A for which the classical path integral 
(30) coincides with the quantum-mechanical path integral of the system, or, in other words, for 
which a theory would possess at the same time deterministic and quantal character. As already 
mentioned, the Grassmannnian ghost variables are responsible for the deterministic nature of the 
partition function. It is obvious that if the ghost sector could somehow be factored out we would 
extend the path integration to all fluctuating paths in q-space. By formally writing 

5 2 A ( SA\ 

Sq k (t)S qi (f) = M M/ ' gTO '^J ' *' , ' m ' n = 1 '-'"' ^ 

we see that the factorization will occur if and only if the (distribution valued) functional Tki{- ■ •) 
is q m independent when evaluated on shell, i.e., Tki(t,t' ,q m ,0) = Fki(t,t'). This is a simple 
consequence of Eq. (20) where the determinant is factorizable if and only if it is q-independent at 
SA/Sq_ = 0. 

In order to provide a correct Feynman weight to every path we must, in addition, identify 

A[q\ = f t2 dt\ m 6 -^, (32) 

Jti d< ?m 

as can be seen from (24) after factoring out the second term. Assuming that L — L(q t , q\) (i.e., a 
scleronomic system) and that the Hessian is regular, the condition (32) shows that Xk = Xk(qi, q'k)- 
In addition, it is obvious on dimensional grounds that [A;] = [q{\. This, in turn, implies that 
Xk = ctkili, where aik is some real (^independent) matrix. To determine the latter we functionally 
expand A in (32) around q^ and compare both sides. The resulting integrability condition reads: 

5 A S 2 A 
{Sji JqW) ^ = ^ qjit) SqWW) ' (33) 

which is evidently compatible with the condition (31). When is diagonalizable we can pass to 
a polar basis and write (32) in more manageable form, namely 

AW = ly^-Mt)^. (34) 

For simplicity, we do not use new symbols for transformed q's. 

To proceed we assume that the kinetic energy is quadratic in q and q. Then Eq.(34) implies 
that Lkin must be liner in q. As such, one can always write (modulo the total derivative) 

ikin = Bij qi(t)qj(i) , (35) 
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with B being an upper triangular matrix. Comparing Lkin on both sides of (34) we arrive at the 
equation 



(a m - l)B im = B mi a m => (B — B T )a = B , (36) 

with no Einstein's summation convention applied here. Because B is upper triangular, the first 
part of Eq.(36) implies that the only eigenvalues of a,j are 1 and 0. Thus, a can be reduced to 
the block form 



a = 



"0 


0" 





1 



(37) 



where 1 is a r x r (r < N) unit matrix. Using the equation (B — B T )a = B we see that B has 
the block structure 



B 






B 2 ' 









(38) 



where B2 is an (N — r) x r matrix. To determine r we use the fact that a is idempotent, i.e., 
a 2 = a. Multiplying (B — B T )a = B by a we find 

Ba = B , B T a = . (39) 

From Ba = B follows that rank(B) = rank(a) = r, whereas B J (I — a) = B T implies that 
rank(£> T ) = rank(l — a). Utilizing the identity rank(£>) = rank(_B T ) we derive r — rank(a) = 
rank(l — a) = (N — r), and thus r = N/2. Thus the condition (34) can be satisfied only for an 
even number N of degrees of freedom. An immediate further consequence of (38) is that we can 
rewrite (35) as 

N/2 

Lkin = ^ B it ( N /2+j) qiqN/2+j ■ (40) 

Denoting ajsr/2+ii QN/2+i arL d ^N/2+i (i — 1, . . . , N/2) as <5j, q~i, and Aj, respectively [hence, A = 
and A — q ], then Eq.(34) reads 

A M = J\ tm ^. (41) 

Here A[<i, q] = A[qi, . . . , Q2n]- The result (41) can be obtained also in a different way. Indeed, in 
Appendix C we show that (34) is a so-called Euler-like functional 



AH - / 

Jti 



« rm\ r-(».W. ■ ■ - r-tftoW. «qg>«M> , . . . , g|eg | , ,42, 



with r(t) being an arbitrary function of % whose variations vanish at the ends 8r{U) = Sr(tf) = 
if all Sqk's have this property. In particular, we may chose r to be any finite power q\! ak (for 
k = 1, . . . , N), in which case 

( - f ^ <?i/? fe " l/afc ) t <^/C /at ) 



^[q]= f'dtq^L 

L J / t fc 1 ai/ak an 

■'ti \ <7i. <7i. 



,0,...,-^ ^ H- (43) 



Assuming, as before, that the kinetic term in L is quadratic in q and q, we arrive at a as in (37), 
and the action (43) reduces again to (41). 

One can incorporate the constraints on on (or Aj) by inserting a corresponding 6- functional into 
the path integral (23). This leads to the most general generating functional with the above-stated 
property: 



Zcm = / X>qX>qX>AX>A 6[X]S[X - q] exp 



= / 2?q2?q exp 



. , &4[q,q] 
i> I dt q — — — 



*2 W 



dty^J k q k 



1 k=i 



/f^ 2 f 
2?q2?q exp i / dtL + dt 



(44) 



An irrelevant normalization factor has been dropped. The Lagrangian L coincides precisely with 
the Lagrangian (3), and describes therefore 'tHooft's deterministic system. Hence within the 
above assumptions there are no other systems with the peculiar property that their full quantum 
properties are classical. Among other things, the latter also indicates that the Koopman-von Neu- 
mann operatorial formulation of classical mechanics [38] when applied to 't Hooft systems must 
agree with its canonically quantized counterpart. 



IV. 'THOOFT'S INFORMATION LOSS AS A FIRST-CLASS PRIMARY CONSTRAINT 

As observed in Section II, the Hamiltonian (1) is not bounded from below, and this is true for 
any function /j. Thus, no deterministic system with dynamical equations qi = /j(q) can describe 
a physically acceptable quantum world. Its Hamiltonian would not be stable and we could build 
a perpetuum mobile. To deal with this problem we will employ 'tHooft's procedure [13]. We 
assume that the system (1) has n conserved, irreducible charges C,, i.e., 

{d,H}=0, t = l,...,n. (45) 

In order to enforce a lower bound upon H , 't Hooft split the Hamiltonian as H = H + — H_ with 
both H + and H- having lower bounds. Then he imposed the condition that if_ should be zero 
on the physically accessible part of phase space, i.e., 

« 0. (46) 

This will make the actual dynamics governed by the reduced Hamiltonian H+ which is bounded 
from below, by definition. 

To ensure that the above splitting is conserved in time one must require that {H-,H} = 
{H + ,H} = 0. The latter is equivalent to the statement that {H + ,H-} = 0. Since the charges 
Ci in (45) form an irreducible set, the Hamiltonians H + and 7J_ must be functions of the charges 
and H: H + = F + (Ck,H) and H- = F_(Ck, H). There is a certain amount of flexibility in finding 
F_ and F + , but for convenience's sake we confine ourselves to the following choice 

[g + Ei°j(*)Ci] 2 [ff-E^iWc.] 2 

+ 4E i Oi(t)C i ' " 4E^(*)Ci ' 1 ' 

where asj(i) arc independent of q and p and will be specified later. The lower bound is then 
achieved by choosing Ej &i(t)Ci to be positive definite. In the following it will also be important 
to select the combination of Cj's in such a way that it depends solely on q (this condition may 
not necessarily be achievable for general / a (q))- Thus, by imposing H- - Owe obtain the weak 
reduced Hamiltonian H w H + w Ej a-i{t)Ci. 

The constraint (46) (resp (47)) can be motivated by dissipation or information loss [14,15,19]. In 
Appendix D we show that the explicit constraint (46) does not generate any new (i.e., secondary) 
constraints when added to the existing constraints (5). In addition, this new set of constraints 
corresponds to 2iV second-class constraints and one first-class constraint (see also Appendix D). It 
is well known in the theory of constrained systems that the existence of first-class constraints signals 



9 



the presence of a gauge freedom in Hamiltonian theory. This is so because the Lagrange multipliers 
affiliated with first-class constraints cannot be fixed from dynamical equations alone [32] . The time 
evolution of observable (physical) quantities, however, cannot be affected by the arbitrariness in 
Lagrange multipliers. To remove this superfluous freedom that is left in the formalism we must 
pick up a gauge, i.e., impose a set of conditions that will eliminate the above redundancy from the 
description. It is easy to see that the number of independent gauge conditions must match the 
number of first-class constraints. Indeed, the requirement on a physical quantity (say /) to have 
a unique time evolution on the constraint submanifold M., i.e., 

m m' 

.1 » {f,H} + J>{/,^} + (48) 

i=l fc=l 

implies that 

{f,<Pi} » 0. (49) 

The constraints p>i and (j)k represent first and second-class constraints, respectively. First-class 
constraints have, by definition, weakly vanishing Poisson's brackets with all other constraints; any 
other constraint that is not first class is second-class. While the Lagrange multipliers Uk can be 
uniquely fixed from the dynamics by consistency conditions (c.f. Appendices A and D) this cannot 
be done for the v^s. In this way (49) represents an obligatory condition for a quantity / to be 
observable. Equation (49) can be considered as a set of m first-order differential equations on the 
constrained surface with the relation {ip i} ipj} w serving as the integrability condition [32,39]. 
Thus, / is uniquely defined by its values on the submanifold of the initial conditions for Eq.(49). 
As a result, the above initial value surface describes the true degrees of freedom. By denoting 
the dimension of the constraint manifold as D we see that the dimension of the submanifold of 
initial conditions must be D — m. We can take this submanifold to be a surface F* specified by 
the equations 

ifii = 0, i = l,...,m, 
4>k = 0, k = 1, . . . ,m! , 

Xi = 0, i = l,...,m. (50) 

The m subsidiary conditions \i are the sought gauge constraints. The functions \i must clearly 
satisfy the condition 

det\\{ X i,Vi}\\ + 0, (51) 

as only in such a case we can determine specific values for the multipliers vt from the dynamical 
equation for \i (this is because the time derivative of any constraint, and hence also \U must be 
zero). Therefore only when the condition (51) is satisfied do the constraints (50) indeed describe 
the surface of the initial conditions. 

The preceding discussion implies that in our case the surface L* is defined by 

<^(q,q,P,p) = 0, x(q,q,P,p) = 0, (52) 
0i(q,q,P,p) = 0, *= 1 2JV. (53) 

The explicit form of p is found in Appendix D where we show that ip w H — J^diCi. Apart from 
condition (51) we shall further restrict our choice of \ to functions satisfying the simultaneous 
equations 

{xAi} - 0, i = l,...,2N. (54) 

Such a choice is always possible (at least in a weak sense) [30] and it will prove crucial in the 
following. 
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In order to proceed further we begin by reexamining Eq.(44). The latter basically states that 



Z C m = J T>qd[q- q c ] exp 



dt q(t)J(t) 



(55) 



We may now formally invert the steps leading to Eq.(8), i.e., we introduce auxiliary momentum 
integrations and go over to the canonical representation of (55). Correspondingly Eq.(55) can be 
recast into 



, 2N 



Zcm = f VpVqDpVq^ |det |{&, <^}||| f[ exp if dt [pq + pq - H] + f dt [qJ + qJ] 

Due to 8- functions in the integration we could substitute 't Hooft's Hamiltonian H for the canonical 
Hamiltonian H. It should be stressed that despite its formal appearance and the phase-space 
disguise, the latter is still the classical partition function a la Gozzi et al. 

To include the constraints (52) into (44) we must be a bit cautious. A naive intuition would 
dictate that the functional 5 functions S[x] and 8[ip] should be inserted into the path-integral 
measure for Zqm- This would be, however, too simplistic as a mere inclusion of S functions into 
Zcm would not guarantee that the physical content of the theory that resides in the generating 
functional Zqm is independent of the choice x- Indeed, utilizing the fact that the generators of 
gauge transformations are the first class constraints [39] we can write that 

$X = £{x,v} + Ctp w e{x,(f}. (56) 

Here e is an infinitesimal quantity. The corresponding gauge generator sip generates the infinites- 
imal canonical transformations 

q^q + <5q, p^p + (5p, 5q = {sip, q} , p = {e^,p}, 

q^q + <5q, p ^ p + <5p , 5q = {eip, q} , p = {e^,p}. (57) 

It follows immediately that the corresponding generating function is 

G(q,q,P,P) =qP + qP + e^ + o(e 2 ). (58) 

The canonical transformations (57) result in changing <p and <pi by 

Sip = Aip, (59) 
6(j>i = e{<pi,p} = B { ip + D i:j <j>j . (60) 

Here A, Bi,C and are some phase-space functions of order e. Note that in our case the gauge 
algebra is Abelian 3 . As a consequence of (59) and (60) we find 

%] ll + Tr^r 1 ^], (61) 
- ll + Tr^)!- 1 ]^^], (62) 

i i 

^det - |l + Tr(£>)| ^dct ||{^,^}|||. (63) 

[here Tr(A) = ^2 t A(t), etc.] In (63) we have used the fact that in the path- integral measure are 
present S[ip] and S[<pi], and so we have dropped on the RHS's of (61)-(63) the vanishing terms. 
The infinitesimal gauge transformations described hitherto clearly show that Zqm is dependent on 
the choice of \ [the term with 1 1 + Tr(^4) | does not get canceled] . To ensure the gauge invariance 



3 If T is any phase-space function then [S e , 8ri\F = S^S^F — 8 v 8 £ J r = e-q {T, {ip, ip}} — 0. 
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we need to factor out the "orbit volume" from the definition of Zcm- This will be achieved by a 
procedure that is akin to the Faddeev-Popov-De Witt trick. We define the functional 



(64) 



with x 9 representing the gauge transformed x- The superscript g in Eq.(64) denotes an element 
of the Abelian gauge group generated by ip. We point out that the functional (64) is manifestly 
gauge invariant since 

(A x v) _1 = Jv 9 5[x a ' 9 ] = J V{g'g)5[x 3 ' 9 ] - (A*)" 1 . (65) 

The second identity holds because of the invariance of the group measure under composition, i.e., 
Vg = V(g'g). Equations (64) and (65) allow us to write "1" as 

1 = A x 5k] J Vg. (66) 
To find an explicit form of A[x] we can apply the infinitesimal gauge transformation (56). Then 
X 9 = X + Z{X,V} + C? => (A^ 1 = J Ve 5[ X + e{ X ,v} + C<p] , 

=> (Axr 1 ^ = IdetlKxMir 1 , (67) 

with the obvious notation det ||{x(i)> = lit {%(*)> Upon insertion of Eq.(66) into 

Zcm we obtain 



Zcm = [ VpVqDpVq |det ||{ X , <p}[\ y/\deb \\{<j>i, <Pj}\\\ 5\x]5[v>]f[5[^ 

rt-2 />t2 

i dt [pq + pq - H] + / dt [qJ + qJ] 



x exp 



(68) 



where the group volume Gy = J T>g has been factored out as desired. The partition function 
(68) is now clearly (locally) independent of the choice of the gauge constraints x- This is because 
under the transformation (59) we have 

det||{xM« - (l + Tr(A))det|{x + <Sx^}|, (69) 

and hence the partition function Zcm as obtained by (68) takes the same form as the untransformed 
one, but with x replaced by x + Because we deal with canonical transformations it is implicit 
in our derivation that the action in the new variables is identical, to within a boundary term, 
with the original action. In path integrals this might be invalidated by the path roughness and 
related ordering problems 4 . For simplicity's sake we shall further assume that the latter are 
absent or harmless. This happens, for instance, when canonical transformations are linear. In 
such cases an infinitesimal change in x does not alter the physical content of the theory present 
in Zcm- This conclusion may generally not be true globally throughout phase space. Global 
gauge invariance, however, is mandatory in our case since we need a global equivalence between 
the partition functions Zcm and Zqm and not mere perturbative correspondence. Thus the 
potentiality of Gribov's copies must be checked in every individual problem separately. 

In passing we may notice that if we arrange the constraints in one set {r] a } — {Xi Pi 4>i} we carL 
write (68) as 



4 In the literature this phenomenon frequently goes under the name of the Edwards-Gulyaev effect [40] . 
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„ 2N+2 

Z C M = / PpPqPpPq y/\det \\{Va,Vb}\\\ 1] %a 
J a=l 



x exp 



dt [pq + pq - H] 



t 2 



dt [qJ + qJ] 



(70) 



By comparison with (7) we retrieve a well known result [39,41], namely, that the set {rj a } of 2N + 2 
constraints can be viewed as a set of second-class constraints. Thus, by fixing a gauge we have 
effectively converted the original system of 2N second-class and one first-class constraints into 
2N + 2 second-class constraints. 

In view of (6) and (54), we can perform a canonical transformation in the full phase space in 
such a way that the new variables are: P\ — \, Qi+i = <p2i, Pi+i = <p2i-i', i = 1,. ..,N. After a 
trivial integration over P a and Qi+i we find that 



ZcM 



-I 



VPVQVQ 1 [6[<p] 



det 



Sip 



SQi 



) eX 4// 



dt 



PQ K 



+ 



r 



dt Qj 



(71) 



where P a and Q a are the remaining canonical variables spanning the (2N — 2)-dimensional 
phase space. To within a time derivative term the new Hamiltonian is done by the prescrip- 
tion K(P, Q, Qi) — H (P, Q, Pi — 0, Qi, Qi+i = 0, P\+i = 0). The sources j are correspondingly 
transformed sources J and J. Utilizing the identity 



6[<p] 



det 


Sip 




SQi 



5[Qi-QI(P,Q)]. 



we can finally write 

ZcM 



I 



VPVQ exp 



*2 



dt 



PQ K* 



rt2 

/ dt Qj 

Jti 



(72) 



(73) 



Here K*(P, Q) = K (P, Q, Qi = Q*(P, Q)). In view of (146) we can alternatively write Z C m as 



Zcm 



I 



VPVQ exp 



dt 



PQ H* + 



dt Qj 



(74) 



where = H + (P, Q, Qi = Ql(P, Q), P a = 0, Qi+i = 0). In passing we may notice that P a and 
Q a are true canonical variables on the submanifold T* of the initial conditions for Eq.(49). Indeed, 
in terms of a non-canonical system of variables {Q} — {ip; \; 4>i\ Q; P} the Poisson bracket of any 
two observable quantities (say / and g) on the constraint manifold M is 



{f,g}\ 



M 



df dg 
d( a dCb 

with {Qj} = {Q;P} and with 



a.b 



df* dg* 
dP t dQ 3 



df* dg* ,_ r . 
2^ Q « 7^7^-' ( 75 ) 



dQ t dQi 



/*(Q,P) = /(^ = 0,x = 0,^ = 0,Q,P), 
5 *(Q,P) - ff (^ = 0,x = 0,<A i = 0,Q,P), 

representing the physical quantities on A4. The latter depend only on the canonical variables Q 
and P which are the independent variables on F*. In deriving (75) we have used the fact that 
various terms are vanishing on account of Eqs.(49) and (54). So, for instance, [{<p, Q} df /dQ]\M = 
0, {(fi,Pj} = 0, {tpi,Qj} — 0, [{XjCj} df/dx}\M — 0, etc. The matrix Ojj stands for the 
(2iV - 2) x (2A^ - 2) symplectic matrix. 

Zqm as defined by (73)-(74) does not generally represent a (classical) deterministic system. This 
is because the constraint <p = explicitly breaks the BRST invariance of Zcm which (as illustrated 
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in Section III) is key in preserving the classical nature of the partition function. Indeed, using the 
relations {x,Pa} = {x>Pa — q a } = 0we immediately obtain 

which implies that 



r xi sr (d X *dtp* d X *d(f*\ _ 

Here the notations X*(q>A) = x(Q:P = A,q = A,p = 0) and <^*(q, A) = ip(q, A, A, 0) were used. 
We also took advantage of the fact that q = A as indicated in Section III. So the generating 
functional (73) (or (74)) can be rewritten as 

Z CM [J = 0} = JvqDXVcVc exp[z5] |det ||{ X W}|| , (78) 

where the integration over the ghost fields was reintroduced for convenience. By reformulating Zqm 
in terms of q, A, c and c we can now easily check the BRST invariance. The BRST transformations 
(25) imply that 

c * 9( P* - -n 

"BRST <p = -g^- eci = -ei.x QBRST <p , 

fe RST V* = set = ~e£x^ BRST <P* ■ (79) 

Here £x aBRST and £x-q BRST represent the Lie derivatives with respect to flows generated by the 
BRST and anti-BRST charges, respectively. Analogous relations hold also for x* ■ Correspond- 
ingly, to the lowest order in e we can write 

<%*] |1 - Tr(^^ ffiBRST )|- 1 J[ X *] , 
|det||{xW}||| - |l-Tr(^ CBRST )||det||{ X *,^}||| . (80) 

The transformations (80) show that the term 5[x*] |det ||{x*, <p*}\ \ in (78) is the BRST invariant 
(as, of course, are both the integration measure and the effective action S). However, because the 
variation i5brst<5[<^*] is not compensated in (78) we have in general, <5brst^cm[J = 0] ^ 0. An 
analogous result applies also to the anti-BRST transformation. 

We should note that the condition i5brst-^cm[J = 0] ^ only indicates that the classical path- 
integral structure is destroyed; it does not, however, ensure that the ensuing Zqm can be recast 
into a form describing a proper quantum-mechanical generating functional. The straightforward 
path-integral representation such as (73) emerges only after the gauge freedom inherent in the 
"information loss" condition ip is properly fixed via the gauge constraint x- Let us finally empha- 
size once more that the partition function (73) (resp. (74)) has arisen as a consequence of the 
application of the classical Dirac-Bcrgmann algorithm for singular systems to the classical path 
integral of Gozzi et al. 



V. EXPLICIT EXAMPLES 



A. Free particle 

Although the preceding construction may seem a bit abstract, its implementation is quite 
straightforward. Let us now illustrate this with two systems. As a warm-up example we start 
with the Hamiltonian 



H = L 3 = xp y - yp x , 
14 



(81) 



which is known to represent the angular momentum with values unbounded from below. Alter- 
natively, (81) can be regarded as describing the mathematical pendulum. This is because the 
corresponding dynamical equation (2) for q is a plane pendulum equation with the pendulum 
constant l/g = 1. The Lagrangian (3) reads 

I = xx + yy + xy - yx . (82) 

It is well-known [42] that the system has two (functionally independent) constants of motion - 
Casimir functions. For (81) they read 

C\ = x 2 + y 2 , C 2 = xp x + yp y . (83) 

The charge C\ corresponds to the conserved radius of the orbit while C 2 is the Noether 
charge of dilatation invariance of the Lagrangian (82) under the transformations (x,y,x,y) 1— » 
(e~ s x,e~ s y,e s x,e s y). As only C\ is p- independent, the functions F + and F_ of this system are 
according to Eq. (47) chosen as: 

(g + aid) 2 (H-a.C,) 2 
F+ 4a 1 C 1 ' F - - 4 ai d ' (84) 

Hence H_ = implies that H + w cii(x 2 +y 2 )- Here a\ is some constant to be specified later. The 
ensuing first-class constraint is 

ip = xp y - yp x - a x x 2 - aiy 2 - p x y + 2a x p x x + p v x + 2aip v y s=s H - a\C\ . (85) 

The gauge condition can then be chosen in the form \ = Py ~ V- Indeed, we easily find that 

{X, f} = Px-x ^ , 

{ X ,<f>i} = 0, i = l,...,4. (86) 

The advantage of our choice of \ is that it will not run into Gribov ambiguities, i.e., the equation 
ip = will have globally unique solution for Qi on T* . This should be contrasted with such choices 
as, e.g., x — Px or x ~ Py, which also satisfy the conditions (86), but lead to two Gribov copies 
each. 

With the above choice of x we ma Y directly write the canonical transformations: 



Pi 


= X 


= Py-y, 


Qi 


= Py 


P2 


= Px 


-X, 


Q2 


= Px 


P3 


= Py 


-y, 


Q 3 


= Py 


p 


= Px 


-x, 


Q 


= Px 



(87) 

It might be checked that the transformation Jacobian is indeed 1. In the new canonical variables 
the Hamiltonian K reads 

K(P,Q,Qi) = H{P,Q,P a = 0,Q 1 ,Q 2 = 0,Q 3 = 0) = -PQ X . (88) 

The functional (5-function (72) has the form 

S[Qi-Qi(P,Q)] = tf[Qi + aiP], (89) 

and hence K*(P,Q) = H+(P,Q) = a\P 2 . Let us now set a\ = l/2mK. After changing variables 
Q(t) to Q(t)/h we obtain not only the correct "quantum- mechanical" path-integral measure 

vQV p K n , (90) 

but also the prefactor l/h in the exponent. So (74) reduces to the quantum partition function for 
a free particle of mass to. As the constant a\ represents the choice of units (or scale factor) for 
C\ we see that the quantum scale h is implemented into the partition function via the choice of 
the "loss of information" constraint. 
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B. Harmonic oscillator 



The system (81) can also be used to obtain the quantized linear harmonic oscillator. This 
is possible by observing that not only C\ = x 2 + y 2 is a constant of motion for (81) but also 
C\ = x 2 + y 2 + c with c being any q and p independent constant. So in particular we can choose 
c = c(q). The functional dependence of c on q cannot be, however, arbitrary. The requirement 
that 't Hooft's constraint should not generate any new (i.e., secondary) constraint represents 
quite severe restriction. Indeed, in order to satisfy Eq.(141) the following condition must hold (c.f. 
Appendix D): 

2N _ 

5>{*,H} = -J2^Pa}{Pa,H} = J>^&^ (91) 

i—0 a,i i,k,a 

which for the system in question is weakly zero only if 

The latter equation has the solution (modulo irrelevant additive constant) c(q) = d 2 (x 2 + y 2 ). 
Here d 2 represents a multiplicative constant. Hence we have that C\ has the general form 

d = x 2 +y 2 +d 2 (x 2 +f). (93) 

It will be further convenient to choose a\ = —l/2d. The resulting first-class constraint then reads 

ip = xp y - yp x + ^x 2 + ^y 2 - ^x 2 - - yp x + xp v - ^xp x - ^yp v + dxp x + dyp y 

«H + ^ Cl . (94) 
If we choose the gauge condition to be 

X = Py + dp x -y, (95) 

it ensures that 

{X, = 2 Ps -2x- 2dp y ^ , 

= 0, i=l,...,4. (96) 

In addition, we shall see that (95) guarantees the unique global solution of the equation <p = for 
Qi on T* (hence it avoids the undesired Gribov ambiguity). 

The canonical transformation discussed in Section IV now takes the form 

Pi = X = Py + dp x -y, Qi = p y , 
Pi = Px - x , Q 2 = p x , 

P 3 = Py-V, Q3 = Py, 

P = p x + dp y -x, Q = Px, (97) 

and the Hamiltonian K reads 

K(P,Q,Q{) = -PQi+dQl-dQ 2 . (98) 
The functional <5-function (72) now has the form 

S[Qi-Q*(P,Q)] = 5 ^-Yd P] - (99) 

This finally implies that the Hamiltonian on the physical space T* has the form K*(P,Q) = 
H+(P, Q) = -(l/4c0P 2 — dQ 2 . By choosing d = —mh/2 and transforming Q i-» Q/K in the path 
integral (73) (resp. (74)) we obtain the quantum partition function for a system described by 
the Hamiltonian: (l/2m)P 2 + (m/2)Q 2 , i.e., the linear harmonic oscillator with a unit frequency. 
This is precisely the result which in the context of the system (81) was originally conjectured 
by 't Hooft in Ref. [14]. Note again that the fundamental scale (suggestively denoted as h) was 
implemented into the theory via the "loss of information" condition. 
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C. Free particle weakly coupled to DufHng's oscillator 



There is no difficulty, in principle, in carrying over our procedure to non-linear dynamical 
systems. As an illustration we will consider here the Rossler system. This is a three-dimensional 
continuous-time chaotic system described by the three autonomous nonlinear equations 

dx 

dt y 

-dt= X + Ay ' 

dz 

— = B + xz -Cz, (100) 

where A, B, and C are adjustable constants. The associated 'tHooft Hamiltonian reads 

H = - Px (y + z) + p y (x + Ay) + Pz (B + xz - Cz) , (101) 

and the Lagrangian (3) has the form 

L = xx + yy + zz + x(y + z) - y(x + Ay) - z(B + xz + Cz) . (102) 

The Rossler system is considered to be the simplest possible chaotic attractor with important 
applications in far-from-equilibrium chemical kinetics [43] . It also frequently serves as a playground 
for studying, e.g., period-doubling bifurcation cycles or Feigenbaum's universality theory. For 
the sake of an explicit analytic solution we will confine ourselves only to the special case when 
A = B = C = 0. With such a choice of parameters the Rossler system can be expressed in a 
scalar form as V = yy + yy — y which ensures its integrability [44] . The latter implies that in this 
regime Rossler's system does not posses chaotic attractors. 

To proceed further, we should realize that because Cj are supposed to be p-independent their 
finding is equivalent to specifying the first integrals of the system (100) (i.e., functions that are 
constant along lines of (x,y,z) satisfying (100)). In other words, the differential equations (100) 
represent a characteristic system for the differential equation {H, C{\ = 0. It is simple to see that 
the first integrals of the above Rossler system are x 2 + y 2 + 2z and ze~ v , hence we can identify 
C\ and C 2 with 

Ci = (x 2 + y 2 + 2z) 2 , C 2 = z 2 e' 2y . (103) 

The previous choice provides indeed positive and irreducible charges. The first class constraint ip 
then reads 

ip = - p x (y + z) + p y x + p z xz - ai(x 2 + y 2 + 2z) 2 - a 2 z 2 e~ 2v 

- p x (y + zz - 4a lX (x 2 + y 2 + 2z)) + py (x + 4a iy (x 2 + y 2 + 2z) - 2a 2 z 2 e - 2y ) 
+ pz(x - zx + 4 ai (x 2 +y 2 + 2z) + 2a 2 ze- 2y ) , 
w H-cud -a 2 C 2 . (104) 

Explicit values of a\ and a 2 will be fixed in the footnote 5. A little algebra shows that the gauge 
condition \ can De selected, for instance, as 

x = Px-y- (105) 

Such a choice satisfies the necessary conditions 

{X,</?} = Py+Pz+x 7^ 0, = 0, i = 1, ... ,6 . (106) 

The above \ a ls° allows us to perform the following linear canonical transformation: 
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Pi 


= x = Px - y, 


Qi 


= Py> 


P2 


— Px % 1 


Q2 


= Px, 


P3 


= Py - y, 


g 3 


= Py, 


Pi 


= Pz - Z, 


Qi 


= Pz, 


Pi 


= (ps/d - z/d)/V2, 


Qi 


= (2d Pz - p s /c + x/c)/y/2 


P2 


= (2c Px - p g /d + z/d)/V2, 


Q~2 


= (x/c - p s /c)/V2. 



(107) 



Here c and d represent arbitrary real constants to be specified later. The transformation (107) 
secures the unique global solution Qi for tp = on T* . To show this it is sufficient to observe that 
[H — aiCi — a 2 C 2 ]\ T , is linear in Q\. Indeed, 

[H - aid - a 2 C 2 ]\ v , = V2c Q X Q 2 - V2c (Q x - Q 2 )Q 2 Pi + d/c {Pi + P 2 )Pi 

- A (Pi) 2 - B Pi(Q 2 ) 2 - C (Q 2 ) 4 , (108) 

with A = 2d 2 (4ai + a 2 ), B = -8\/2aidc 2 and C = 4aic 4 . As a result 

K*(P,Q)=H* + (P,Q) = A(Pif+BP 1 (Q 2 f+C (Q 2 ) 4 . (109) 

Inserting this into (73) (resp. (74)) and integrating over P\ and P 2 we obtain the following chain 
of identities: 



Zc 



M 



/ 



VPVQ, exp\i dt [PQ-A {Pif-B Pi{Q 2 f-C (Q 2 ) 4 + Qj 



VQ{DQ 2 5[Q 2 ] exp 



«2 



dt 



lim / VQ{DQ 2 exp I i / dt 
J [Ju 

x exp j^ij dt 



(Qi-B(Q 2 ) 2 ) 2 -C(Q 2 ) 4 + Qi 



4A 



c {Q2T + Qi 



(110) 



As an explanatory step we should mention that the formal measure in the second equality of (110) 
has the explicit time-sliced form 



(111) 



while in the third equality the shorthand notation VQ{DQ 2 stands for 

vqivq 2 « n d 4^) 

. \v4:irieA V^niaeJ 



(112) 



The symbol e represents the infinitesimal width of the time slicing. During our derivation we have 
used the Fresnel integral 



j — 1 



dx e 



+ixi _ l*_ p i(f/a-ir)/i _ l*_ if/da) 



a > 0, 



and the ensuing representation of the Dirac (5-function: 



lim 

a^o + V Aiira 



(113) 



(114) 



In the following we perform the scale transformation Q 2 /yfa 1— > \[2~m 2 Q 2 and set A = l/(2mi), 
B = l/(^/mim 2 ) and C — l/m 2 . 5 The resulting partition function then reads 



5 This choice is equivalent to the solution: 
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= g lim J VQ{DQ 2 expjijT'dt^ (4) 2 + ^ (4) 2 ] } 



r*2 

x cxp I i I dt 



(h((hr - ^f- (Q 2 ) 4 + Qj 



(115) 



where we have set g = 2\/2a. The system thus obtained describes a pure anharmonic (Dufhng's) 
oscillator (Q 2 oscillator) weakly coupled through the Rayleigh interaction with a free particle (Qi 
particle). Alternatively, when mi = m 2 — m we can interpret the Lagrangian in (115) as a planar 
system describing a particle of mass m in a quartic scalar potential e$(Q) = mg 2 /4 (Q2) 4 and 
a vector potential eA = (gm-\/l/2 (<32) 2 ,0) (i.e., in the linear magnetic field B 3 = e 3i jdiAj = 
-gmV2 Q 2 /e). 

It is preferable to set mi 1— > mift and m 2 1— » m 2 /h. The latter corresponds to the scale factors 
C12 = l/(2mift) and a\ — l/(8mift). After rescaling Qi(t) 1— » Q\(t)/h the partition function (115) 
boils down to the usual quantum-mechanical partition function with the path-integral measure 

™ - IT ( dQl{ti) dQ2{ti) ) . (us) 

. y y / 2irieh/mi y / 2itieh/m2 I 

and with l/7i in the exponent. Hence, just as found in the previous two cases, the choice of 
't Hooft's condition ensures that the Planck constant enters the partition function (115) in a 
correct quantum-mechanical manner. In turn, h enters only via the scale factors ai and a 2 (the 
factors d and c are h independent) and hence it represents a natural scale on which the "loss 
of information" condition operates. In other words, whenever one would be able to "measure" 
or determine from "first principles" the "loss of information" condition one could, in principle, 
determine the value of the fundamental quantum scale h. 

As a final note we mention that the 't Hooft quantization procedure can be straightforwardly 
extended to other non-linear systems and particularly to systems possessing chaotic behavior (e.g., 
strange attractors). In general cases this might be, however, hindered by our inability to find the 
corresponding first integrals (and hence Cj's) in the analytic form. It is interesting to notice that 
machinery outlined above allows to find the emergent quantistic system for the configuration-space 
strange attractors. This is because in 't Hooft's "quantization" one only needs the dynamical 
equations in the configuration space. The latter should be contrasted with the Hamiltonian (or 
symplectic) systems where strange attractors cannot exist in the phase-space on account of the 
Liouville theorem [45]. 



VI. CONCLUSIONS AND OUTLOOK 



In this paper we have attempted to substantiate the recent proposal of G.'t Hooft in which 
quantum theory as viewed as not a complete final theory, but is in fact an emergent phenomenon 
arising from a deeper level of dynamics. The underlying dynamics are taken to be classical 
mechanics with singular Lagrangians supplied with an appropriate information loss condition. 
With plausible assumptions about the actual nature of the constraint dynamics, quantum theory 



ai 1 1 

fl i = —r i "= — , c — +— — 
4 2V2a 2 mi 

Without loss of generality we can set d — 1/2, then: 

1 1 1 o3/4 4 / m l 

0-2 = 75 , Ol = , C = ±2 ' . 4 / . 

2mi 8mi y rri2 
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is shown to emerge when the classical Dirac-Bergmann algorithm for constrained dynamics is 
applied to the classical path integral of Gozzi et al.. 

There are essentially two different tactics for implementing the classical path integrals in 
't Hooft's quantization scenario. The first is to apply the configuration-space formulation [26]. 
This is suited to situations when 't Hooft's systems are phrased through the Lagrangian descrip- 
tion. The alternative approach is to start with the phase-space version [27]. The latter provides 
a natural framework when the Hamiltonian formulation is of interest or where the language of 
symplectic geometry is preferred. It should be, however, stressed that it is not merely a matter 
of a computational convenience which method is actually employed. In fact, both approaches are 
mathematically and conceptually very different (as they are also in conventional quantum mechan- 
ics [11,46]). Besides, the methodology for handling singular systems is distinct in Lagrangian and 
Hamiltonian formulations (c.f. Rcfs. [39,41] and citations therein). In passing, we should mention 
that the currently popular Hamilton- Jacobi [47] and Legendre-Ostrogradskii [48] approaches for a 
treatment of constrained systems, though highly convenient in certain cases (e.g., in higher-order 
Lagrangian systems), have not found as yet any particular utility in the present context. 

Throughout this paper we have considered only the configuration-space formulation of classi- 
cal path integrals. (Incidcntly, the phase-space path integral which appears in Section IV (after 
Eq.(55)) is not the phase-space path integral a la Gozzi, Reuter and Thacker [27] but rather 
Gozzi's configuration-path [26] integral with extra degrees of freedom.) By choosing to work 
within such a framework we have been able to render a number of formal steps more tractable 
(e.g., BRST analysis is reputed to be simpler in the configuration space, uniqueness proof for 
't Hooft systems is easy and transparent in the Lagrange description, etc.). The key advantage, 
however, lies in two observations. First, the position-space path integral of Gozzi et al. provides 
a conceptually clean starting point in view of the fact that it represents the classical limit of 
both the stochastic-quantization path integral and the closed-time-path integral for the transition 
probability of systems coupled to a heat bath. Such a connection is by no means obvious in 
the canonical path-integral representation as both the Parisi-Wu stochastic quantization and the 
Fcynman- Vernon formalism (with ensuing closcd-time-path integral) are intrinsically formulated 
in the configuration space. Second, according to 't Hooft's conjecture the "loss of information" 
condition should operate in the position space where it is supposed to eliminate some of the 
transient trajectories leaving behind only stable (or near to stable) orbits [14]. Hence working in 
configuration space may allow one to probe the plausibility of 't Hooft's conjecture. The price that 
has been paid for this choice is that the configuration space must have been doubled. This is an 
unavoidable step whenever one wishes to obtain first-order autonomous dynamical equations di- 
rectly from the Lagrange formulation (a fact well known in the theory of dissipative systems [49]). 
Our analysis in Appendix BII suggests, that the auxiliary coordinates q~i may be related to rel- 
ative coordinates on the backward-forward time path in the Feynman- Vernon approach. (Such 
coordinates also go under the names fast variables [50] or quantum noise variables [51].) On the 
formal side, the auxiliary variables q- L are nothing but Gozzi's Lagrange multipliers A, (in our case 
denoted as Aj). 

In order to incorporate the "loss of information" into our scheme, we have introduced in Section 
IV an auxiliary momentum integration to go over to the canonical representation. Such a step, 
though formal, allowed us to treat our constrained system via the standard Dirac-Bergmann 
procedure. It should be admitted that such a choice is by no means unique - e.g., methodologies 
for treatment of classical constrained systems in configuration space do exist [39,41]. The decision 
to apply the Dirac-Bergmann algorithm was mainly motivated by its conceptual simplicity and 
direct applicability to path integrals. On the other hand, we do not expect that the presented 
results should undergo any substantial changes when some another scheme would be utilized. It 
should be further emphasized that while we have established the mathematical link (Eqs.(52) and 
(146)) between the "loss of information" condition and first-class constraints, it is not yet clear 
if this connection has more direct physical interpretation (although various proposals exist in the 
literature [24,14,19]). Such an understanding would not only help to develop this approach for 
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more complicated physical situations but also affiliation in a systematic fashion of a quantum 
system to an underlying classical dynamics. Work along those lines is currently in progress. 

To illustrate the presented ideas we have considered two simple systems; the planar pendulum 
and the Rossler system. In the pendulum case we have taken advantage of free choice of an additive 
constant in the charge C\. This in turn, allowed us to imposed 't Hooft's constraints in two distinct 
ways. In the case of Rossler's system two p- independent, irreducible charges C\ and C 2 exist. For 
dcfiniteness sake we have constructed in the latter case the "loss of information" condition with 
the additive constant set to zero. With this we were able to convert the corresponding classical 
path integrals into path integrals describing a quantized free particle, a harmonic oscillator, and 
a free particle weakly coupled to Duffing's oscillator. As a byproduct we could observe that 
our prescription provides a surprisingly rigid structure with rather tight maneuvering space for 
the emergent quantum dynamics. Indeed, when the classical dynamics is fixed, the 't Hooft 
condition is formulated via linear combination of charges Ci which correspond to the first integrals 
of the autonomous dynamical equations for q, i.e., Eq.(2). Due to the explicit form of 't Hooft's 
Hamiltonian the constraint is of the first class and so we must remove the redundancy in the 
description by imposing the gauge condition \. By requiring that the consistency conditions (51) 
and (54) are fulfilled, that the choice of \ does not induce Gribov ambiguity, and that the canonical 
transformations defined in Sec. IV are linear, we substantially narrowed down the class of possible 
emergent quantum systems. Note also, that when we start with the TV-dimensional classical system 
(q variables), the emergent quantum dynamics has N — 1 dimensions (Q variables). Indeed, by 
introducing the auxiliary degrees of freedom q we obtain 47V-dimensional phase space which is 
constrained by 2N + 2 conditions (<fii, cp and x), which leaves behind (2N — 2)-dimcnsional phase 
space Q, P. This disparity between the dimensionality of the classical and emergent quantum 
systems vindicates in part the terminology "information loss" used throughout the text. 

An important conclusion of this work is that 't Hooft's quantization proposal seems to provide 
a tenable scenario which allows for deriving certain quantum systems from classical physics. It 
should be stressed that although wc assumed throughout that the deeper level dynamics is the 
classical (Lagrangian or Hamiltonian) one, there is in principle no fundamental reason that would 
preclude starting with more exotic premises. In particular, our conceptual reasoning would go 
unchanged if we had begun with Lagrangians operating over coordinate superspaces (pseudoclas- 
sical mechanics [52]) or with the currently much discussed discrete classical mechanics (i.e., having 
foam-, fractal-, or crystal- like configuration space) [53], etc. . The only prerequisite for such ap- 
proaches is the possibility of formulating a corresponding variant of Gozzi's path integral, and a 
method for implementing the "loss of information" constraint in such integrals. 

There are many interesting applications of the above method. Applications to chaotic dynam- 
ical systems especially seem quite pertinent. After all, central to our reasoning is a (doubled) 
set of real first-order dynamical equations 6 which, under favorable conditions, may by associated 
with a chaotic dynamics in the configuration space. Wc should emphasize that the reader should 
not confuse the above with the extensively studied but unrelated notion of chaos in Hamiltonian 
systems - we do not deal here with dynamical equations on symplectic manifolds. This is impor- 
tant, as Hamiltonian systems forbid per se the existence of attractive orbits which are otherwise 
key in 't Hooft's proposal. In this respect our approach is parallel with some more conventional 
approaches. Indeed, a direct "quantization" of the equations of motion - originally proposed by 
Feynman [54] - is one of the techniques for tackling quantization of dissipative systems [55,56]. 
In field theories this line of reasoning was recently progressed by Biro, Muller, and Matinyan [19] 
who demonstrated that quantum gauge field theories can emerge in the infrared limit of a higher- 
dimensional classical (non-Abelian) gauge field theory, known to have chaotic behavior [57]. 



Non-trivial are only the equations over actual configuration space. The dynamical equations for the 
auxiliary variables q\ are linear and hence they are not relevant in this connection. 
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We finally wish to comment on two more points. First, in cases where one strives for an 
explicit reparametrization invariance (or general covariance) of the emergent quantum system 
the presented framework is not very suitable. The absence of explicit covariance in both Dirac- 
Bcrgmann and Fadeev-Senjanovic algorithms makes the actual analysis very cumbersome or even 
impossible. In fact, expressions (68) and (70) are evidently not generally covariant due to the 
presence of time-independent constraints in the measure. Although generalizations that include 
covariant constraints do exist [33,58,59] they result in gauge fixing conditions which depend not 
only on the canonical variables but also on the Lagrange multipliers (or explicit time) . Such gauge 
constraints are, however, incompatible with our Poisson bracket analysis used in Section IV, and 
Appendixes A and D. Hence, if the emergent quantum system is supposed to be reparametrization 
invariant (e.g., relativistic particle, canonical gravity, rclativistic string, etc.) a new framework 
for the path-integral implementation of 't Hooft's scheme must be sought. Second, the formalism 
of functional integrals is sometimes deceptive when taken too literally. The latter is the case, 
for instance, when gauge conditions are imposed and/or canonical transformations performed. 
The difficulty involved is known as the Edwards-Gulyaev effect [11,40,46] and it resides in the 
exact nature of the limiting sequence of the finite dimensional integrals which constitute the path 
integral. As a result the classical canonical transformation does not leave, in general, the measure 
of the path integral Liouville invariant but, instead induces an anomaly [46,60]. Thus, for our 
construction to be meaningful it should be shown that the canonical transformations in Section IV 
are unaffected by the Edwards-Gulyaev effect. Fortunately, in cases when the generating function 
is at most quadratic (making canonical transformations linear) and not explicitly time dependent, 
it can be shown [29,60,61] that the anomaly is absent. It was precisely for this reason that more 
general transformations were not considered in the present paper. Clearly, both mentioned points 
are of key importance for further development of our procedure and, due to their delicate nature, 
they deserve a separate discussion. 

Let us end with the remark that the notorious problem with operator ordering known from 
canonical approaches has an elegant solution in path integrals. The ordering is there naturally 
generated by the necessary physical requirement that path integrals must be invariant under 
coordinate transformations [65]. 
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APPENDIX A 

In this appendix we show that the system (1) has no secondary constraints. In contract to the 
primary constraints which are a consequence of the non-invertibility of the velocities in terms of 
the p's and <?'s, secondary constraints result from the equations of motion. To show their absence 
in 't Hooft's system we start with the observation that the time derivative of any function /(q, p) 
is given by [39] 

/ » {f,H} + (117) 

Here u a are the Lagrange multipliers to be determined by the consistency conditions 

« 4>i « {&,H} + (118) 

The latter is nothing but the statement that constraints (as functions of q and p) must hold at 
any time. If all could not be determined from the consistency condition (118) then we would 
have the so-called secondary constraints. In our case we have 
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dH 

dq a 



96 0, {<j>%, H} = -/ n (q) 96 0, = Sab 



(119) 



Using the fact that {<pi,H} 96 and det \{<f>i, 4>j}\ = 1, the inhomogeneous system of linear 
equations (118) can be uniquely resolved with respect to u- 7 , thus implying the absence of secondary 
constraints. 



APPENDIX B 



BI 



We show here that Gozzi's configuration-space path integral results from the "classical" limit of 
the stochastic-quantization partition function, i.e., the limit where the width of a noise distribution 
tends to zero. For this purpose we start with the form of the partition function for stochastic 
quantization as written down by Zinn- Justin [34,62]: 



Z SC (J) = J VqVcVcVX exp |— <S[q, c, c, A] + J J(x)q 



(x)dx 



(120) 



where 



«,(A) + / X(x) (°B& + J± 1 <U- 



\ dr 



5q{x) 



and 



- [ dxdx' c a (x) f-^-S ab 6(x - x') + - / f ) c b (x') 
J \dr 5q a (x)6q b (x') J 

exp[w(X)] = J XVcxp |— a{v) + J dx\(x)v(x)^ , 



(121) 



(122) 



with T>v cxp(— a(y)) being the functional measure of noise. Here x = (i, t) and dx = dtdr where 
r is the Parisi-Wu fictitious time. The dynamical equation for q(x) is described by the Langevin 
equation 



dq{x) SA[q] 



dr 



5q 



(123) 



q=q(x) 



with the initial condition q(t, 0) = q(i). For Gaussian noise of variance 2ft, the noise measure is 
^e X p(-^)) = n^exp(-^ J dx V \x)^ , (124) 

and (120) takes the form 

Z SC (J) = [vqDus(^ + S -^--u)det ^-5 ab S(x - x') + - .**f . ,. 

J \dr Sq ) Ut 5q a {x)dq b (x') 

x cxp|— o(v) + j 3{x)q{x)dx^ 

(125) 



= J VqDv 5 q-q M cxp j-cr(z/) + J 3(x)q{x)dx^ . 
where <5[/(q)] = FJt T S(f(q(t, r))) and qM(x) is a solution of (123). Using the representation. 



S(x) 



lim — == 

2vVft 



-x 2 /(4h) 



(126) 
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we get in the limit of zero distribution width (i.e., h — > 0+) that 

Zsc{J,K) -> I T>q 5\q- q [0] j exp { [ J{x)q(x)dx 



(127) 



Choosing a special source J (a;) = J(t)5(r) we can sum in the path integral solely over configurations 
with q(t, 0) = q(t) as other configurations will contribute only to an overall normalization constant. 
Inasmuch we finally obtain 



lim Z sc {3,h) = Z CM {3). 

h— >o+ 



(128) 



BII 



In this part of the appendix we show that Gozzi's configuration-space partition function (19) 
results from the "classical" limit of the closed-time path integral for the transition probability of 
a system coupled to a thermal reservoir at some temperature T. By the classical limit we mean 
the high temperature and weak heat bath coupling limit. 

The path-integral treatment of systems that are linearly coupled to a thermal bath of harmonic 
oscillators was first considered by Feynman and Vernon [63] . For our purpose it will be particularly 
convenient to utilize the so called Ohmic limit version, as discussed in Refs. [11,64]: 

Z FV [J+,J-] = J Vq+Vq_ exp L4[q+] - .4[q_]] + J dt [J + (i)q+(i) - J_(t)q_(t)] J 

x exp {~^ / dt [^)-q-W][q+W + q-Wp} 

x exp {~WpJ dt J dt ' [q+(*)-q-(*)]^(*'*')[q+(*')-q-(0]} • (129) 

Here the paths q+(£) and q~(t) are associated with the forward and backward movement of 
the particles in time. The super-script R indicates a negative shift in the time argument of the 
velocities with respect to positions. The latter ensures the causality of the friction forces [64]. In 
addition, m represents the particle mass (for simplicity we assume here that all system particles 
have the same mass), (3 — 1/T, and 7 is the friction constant (or thermal reservoir coupling). 
The function K(t,t') is the bath correlation function. As argued in [11,64], at high temperatures 
K(t, t') rts 5(t — t'). Introducing the new set of variables q = [q + + q_]/2 and q = [q + — q_] (i.e., 
the center-of-mass and fast coordinates) we can in the high-temperature case recast (129) into 

Z FV [J,J] = J VqDq cxp|^[^l[q + q/2]-^[q-q/2]] + J dt [J(t)q(t) - J(t)q(t)] J 

X exp j-i^ J dt q(t) [q(t)f - ^1 J dt q 2 (t)| . (130) 

Here the self-explanatory notation J = [J + — J_] and J = — [J+ + J-]/2 was used. Let us now 
define lo — 2mj/(3, integrate over q, and go to the classical limit 7 — > 0. Then we obtain the 
following chain of equations: 
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lim Z FV [J,J] 

7^0 



lim / VqDq exp < — dt q(i) 
7^0 J [hj 



' )A m 1 [q{t)f + iM{t) 



5q(t) 



x exp dt J(t)q(t) 



= lim / Vq exp < / dt 



i|-m 7 [q(t)] fl + l fiJ(t) 



lim yPqJ[q] exp --/dt 



n 2 



= /pq* 
= fvqS 



<5q 



q- q 



[J] 



<y q (t) 

J[q] expjy" dt J(t)q(t)| 
exp i / dt J(t)q(t) > • 



mjq(t) + ihJ(t) 



^/*q a (*) 



dt J(t)q(t) 



dt J(t)q(t) 



(131) 



The Jacobian ,7[q] results from transition to the "unretarded" velocities and its explicit form 
reads [64]: 



J[q\ = det 



9 r r /, ^ 



Coordinates q^ J ' are solutions of the equation of the motion: 

5A[q] 



Sq(t) 



-iftj(t) . 



In the limit 7 — > 0, we find again the Gozzi et al. partition function 

limZpvlJjO] = lim lim Zpv [J j J] = ^cm[J] 

7^0 ' h—>o~/^o 



(132) 



(133) 



(134) 



APPENDIX C 

In this appendix we prove that (34) is a special case of the Euler-like functionals (42). Let 
us first show that (34) can be replaced by an action of the form (42). Indeed, because of the 
homogeneity of (34), we can immediatley replace it by 



« = e/* *^ t)qi ® wfi«( t ) = S dt r{t) 

Since this is true for any r(t), we see that 



6A[r a * qi 
8r(t) 



1 



(135) 



(136) 



This simply expresses the fact that the functional A[r ai qi] is linear in r(t). The right-hand side 
of (135) has then precisely the Euler form (42). 



The reverse direction is proved in the following way: We first recast (42) in the general form 
dtr(f)L(q(t),q(t)) = J dt L(r a * (t) qi (t), d(r a > (t) qi (t))/dt) . (137) 
Applying the variation J dt 8/5r(t) to (137) we obtain 
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This relation must hold for all r(t), and hence by choosing r(t) = 1 we arrive at the required result 

A[q\ = fdt ^^(t)^. (139) 



APPENDIX D 



Here we prove the fact that inclusion of the subsidiary constraint (46) in the primary constraints 
(5) does not produce any secondary constraints. The secondary constraints result from the con- 
sistency conditions (118) or, in other words, when existent constraints are incompatible with the 
equation of motion. 

We first observe that the condition H- w can be equivalently represented by the condition 
(H — J2i a iCi) = fa ~ 0. If we now add the subsidiary constraint fa to the remaining 2N 
constraints fa and again require that the constraints fa remain (weakly) zero at all times we have 







fa w {fa,H} + u j {fa,fa}, i,j = 0,1..., 2N. 



(140) 



Since there is an odd number of constraints and because {fa, fa;} is an antisymmetric matrix we 
have that det || {fa, fa} || = 0. From the analysis in Appendix A it is clear that the rank of the 
matrix {fa, fa} is 2N and hence it has one null-eigenvector, say e. Inasmuch, Eq.(140) implies 
the constraint 



2N 

~Y^ei{fa,H} 



(141) 



i=0 



If the latter would represent a new non-trivial constraint (i.e., constraint that cannot be written 
as a linear combination of constraints fa) we would need to include such a new constraint (the 
so called secondary constraint) into the list of existent constraints and go again through the 
consistency condition (140). Fortunately, the condition (141) is automatically fulfilled and hence 
it does not constitute any new constraint. Indeed, be choosing 



/ 1 \ 

{fa,fa 2 } 

{fa^\} 

{fa\,fa} 

{fa, fa?} 
\{K,M J 



( 1 \ 

/a(q) 

dcf>o 
dq a 

A(q) 

d<t>a 



/iv(q) 



(142) 



and using {fa, H} — together with (119) we obtain 



(143) 



i=l 



As the latter is zero (even strongly) there is no new constraint condition generated by an inclusion 
of fa in the original set of (primary) constraints. Note, that the key in obtaining (143) was the 
fact that Ci's are p-independent constants of motion. 



The rank of {fa, fa;} being 2N means that there is one relation 



2N 



^2ei{fa,fa} w 0. 
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(144) 



i=0 



Any linear combination of the constraints fa is again a constraint. So, particularly if we define 
ip = J2i e i4>i we obtain that <p has weakly vanishing Poisson brackets with all constraints, i.e., 

« 0, i = l,...,2JV. (145) 

Thus, according to Dirac's classification (see e.g., Ref. [32]) is a first class constraint. The 
remaining 2N constraints (which do not have vanishing Poisson brackets with all other constraints) 
are of the second class. Note particularly that the explicit form for <p reads 

v = y>& = (if-v^)-^-,^, (we) 

2—0 1=1 a— 1 

which is clearly weakly identical to H — ^ o^C;. Observe that it is H and not H that is present 
in (146). 
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